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ABSTRACT 

The C/O ratio is predicted to regulate the atmospheric chemistry in hot 
Jupiters. Recent observations suggest that some exo-planets, e.g. Wasp 12- 
b, have atmospheric C/O ratios substantially different from the solar value of 
0.54. In this paper we present a mechanism that can produce such atmospheric 
deviations from the stellar C/O ratio. In protoplanetary disks, different snowlines 
of oxygen- and carbon-rich ices, especially water and carbon monoxide, will result 
in systematic variations in the C/O ratio both in the gas and in the condensed 
phase. In particular, between the II2O and CO snowlines most oxygen is present 
in icy grains - the building blocks of planetary cores in the core accretion model - 
while most carbon remains in the gas-phase. This region is coincidental with the 
giant-planet forming zone for a range of observed protoplanetary disks. Based on 
standard core accretion models of planet formation, gas giants that sweep up most 
of their atmospheres from disk gas outside of the water snowline will have C/0~1, 
while atmospheres significantly contaminated by evaporating planetesimals will 
have stellar or sub-stellar C/O when formed at the same disk radius. The overall 
metallicity will also depend on the atmosphere formation mechanism, and exo- 
planetary atmospheric compositions may therefore provide constraints on where 
and how a specific planet formed. 

Subject headings: astrochemistry — circumstellar matter — planetary systems 
— molecular processes — planets and satellites: atmospheres — planet-disk 
interactions 
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Introduction 



Giant planets are commonly detected around main-sequence stars. A handful of gas 
giants orbiting close to their host stars have been investigated spectroscopically and with 
narrow-band photometry, resulting in constraints on their atmospheric compositions. The 



atmospheres o 



(iFortnev et al. 



2006 



thes e 'hot Jupiters' are predicted to be abundant in H2O, CH4 and CO 



Tinetti et al. 



2010), and observations qualitat i vely confirm these predictions (IRichardson et al 



20071 : iGrillmair et al.l l2008l : ISwain et al.l |2010[ ) . The relative abundances 



of molecule s can potentially be used to characterize atmospheric temperature and pressure 
structures (IFortney et al.l I2OIOI ) as well as to enhance our understanding of atmospheric 
chemistry under extreme conditions. However, chemical predictions are sensitive to devia- 
tions from solar C/0 ratios; for a typical hot Jupiter atmosphere at 1100 K and 0.01 bar, 
an increase in the C/0 ratio from the solar 0.54 to unity may res ult in CH4, C2H2 and CN 
abundance increases by three orders of magnitude (lLoddersll2009l ). 



These extreme effects have recently been used to estimate the C/0 ratio in WASP 12-b. 
Based on comparison between a n arr ay model predictions and six narrow band photome- 
try points, iMadhusudhan et al.l (120 111 ) report an atmospheric C/0~1, despite a solar C/0 
ratio in the host star and an overall planetary metallicity (C/H) that is at most twice the 
stellar value, but more likely sub-stellar (2 x 10^^ — 1 x 10^^ in WASP 12-b versus the 
super-solar metallicity of the host star corresponding to C/H~ 6 x 10~^; 



Madhusudhan et al.ll201ll ). Based on analyzing the observations. IMadhusudhan et al. 



Hebb et a. 



2009 



fl201lh 



further argue that the spectral features originate in the deep layers of the atmosphere, that 
day-night energy redistribution is efficient and that a thermal inversion layer deeper than 
0.01 bar can be ruled out. If the low C/H and high C/0 are confirmed, and the observations 
do trace the bulk of the atmosphere, this planet's composition cannot be explained by previ- 
ously suggested mechanisms to enhance the planetary C/0 ratio. These processes — pile- up 
of carbon-rich grains in the disk, and co ld trapping of H2O (and thus a large fraction of the 



oxygen budge t) outside of the snowline (IStevenson fc LunindllQSSl : iKuchner fc Seagerll2005 
Loddersll2009l ) — cannot account for an enhanced C/0 ratio without also increasing C/H. 



We instead investigate the effects of the different snowlines of a few common volatile 
species as a mechanism for modifying the planet atmosphere C/0 ratio compared to stellar 
values. In this letter we develop a simple, static model to demonstrate the possible effects of 
this mechanism (§2). We discuss possible extensions of the model, especially the inclusion of 
dynamics, in §3. In §4 we explore predictions of planetary atmospheric compositions based 
on the location of the forming planet with respect to different snowlines and the mechanism 
through which the planet atmosphere forms. 
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2. Model prescription and results 



In the core accretion scenario, planet formation begins witli the coagulation of icy grains 
in the protoplanetary disk. Grains agglomerate into planetesimals (km sized bodies), which 
in turn combine t o form the planet c ore. When sufficient ly massive, the core accretes a 
gaseous envelope (jPoUack et al.l llQQGt iHubickyj et al.l |2005| ) . Because the main molecular 
sources of C and O have different condensation temperatures (Table 1), the relative amounts 
of C and O differ in grains and in gas at some disk radii. Atmospheric C/0 need not, 
therefore, reflect the average C/0 of the disk. Instead, the C/0 ratio of a planetary envelope 
depends on where in the disk the planet forms, on the fractions of the atmosphere that are 
accreted from gas and from icy bodies, and on the importance of core dredging, i.e. how 
isolated the atmosphere is from the core. In the simplest case, the core and atmosphere 
are completely isolated from each other, and the atmosphere is built up purely from gas. 
We therefore begin with only considering gas accretion, and then show how the expected 
atmosphere composition is modified by adding planetesimal accretion. 

Once a core is massive enough to begin runaway accretion of a gas envelope, this ac- 
cretio n most likely happen s faste r than the planet can migrate due to interactions with the 
disk. iD'Angelo fc Lubowl (120081 ) estimate that a planet migrates inward by <20% of its 
semi-major axis during runaway growth. We therefore assume that the planetary envelope is 
accreted between the same set of snowlines where accretion started. As a first step, we fur- 
ther assume that grains contributing to the atmosphere come from the same location as the 
gas (which need not be the case) and that gas and grain compositions are constant between 
each set of snowlines. Finally we assume that the snowlines are static, which is justified by 
the long timescaies at which disk midplane temperatures change in disks older than 10^ years 
(when gas giants are proposed to form) compared to the 10^ y ear timescaies of runaway gas- 
accretion ( jLissauer et al.l 120091 : iDodson- Robinson et al.ll2009l ). Specifically, the temperature 
structure is set b y viscous dissipation in the inner disk and irradiation by the central star 
in the outer disk (ID'Alessio et al.lll998l ). and both accretion and stell ar luminosity decay on 



10^ year timesca ies at the time of planetary envelope accretion (e.g. iHartmann et al.lll998 
Siess et al.ll2000l ). We return to these considerations in §3. 



We estimate the total abundances (grain + gas) of the major O- and C-cont aining species 



in typ ical disks from a combination of ice observations of a protoplanetary disk (jPontoppidan 



20061 ). and derived grain compositions in the dense interstellar medium (ISM) (Table 1). The 



fWhittet 


2010 


). 


1 grains 


Draine 



latter carbon sources are unknown, and a high evaporation temperature is adopted to pre- 
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vent this unknown carbon component from influencing the model outcome; if any of this 
carbon is present in more volatile forms, it will enhance the gas-phase C/0 ratio further. 
The sublimation temperature for silicate grains are set to 1500 K. For all other molecules, 
we calculate the d ensity -dependent sublimation temperatures following the prescription of 
HoUenbach et all toO^ using b inding energi es of H^O, CO? and CO of 5800 K, 2000 K 
and 850 K (ICoUings et all |2004J : iFraser et al.l I2OOII : lAikawa et all 119961 ) . A complication is 
the observed ease with which H2O can trap other molecules in its ice m atrix. It is however 
difficult to trap more than 5-10% of the total CO abundance in H2O ice (IFayolle et al.ll201ll ) 
and we therefore ignore this process. 

The radii of different snowlines are set by the disk temperature profile. Co nsistent with 
the temperature profile derived from the comp ositions of solar system bodies ( Lewis! Il974l ) 
and with observations of protoplanetary disks (lAndrews &: Williams! !2005! . !2007! ) we adopt 
a power law profile. 



where Tq is the temperature at 1 AU and q is the power law index. In a large sample o f 



protoplanetary disks, the average Tq is 200 K, and q = 0.62 ( lAndrews &: Williams! !2007l ). 
Figure 1 displays the C/0 in the gas and in grains in the disk midplane as a function 
of distance from the young star for this average disk profile. Between the H2O and CO 
snowlines, the gas-phase C/0 ratio increases as 0-rich ices condense, with the maximum 
C/0~1 reached between the CO2 and CO sublimation lines at 10-40 AU. In the case of 
completely isolated core and atmosphere accretion, the atmospheric C/0 ratios will refiect 
the gas phase abundances, resulting in C enrichments beyond the H2O snowline. 

The size and position of the disk region where the C/0 ratio in the gas reaches unity 
depend on the disk temperature profile. A more luminous star will heat the disk further, 
pushing the various snowhnes outwards, while the steepness of the disk temperature profile 
determines the spaci ng of the different snowlines . Figure 2 compares protoplanetary disk 



thermal profiles from lAndrews &: Williams! ( 12005! ). which sample stars with a range in spec- 



tral types, with the 'typical' disk profile from Fig. 1. In all cases, the gas-phase C/0 ratio 
is enhanced in regions associated with gas-giant formation, i.e. a few to a few tens of AU. 
Formation of C-rich atmospheres from oxygen depleted gas accretion can therefore operate 
in most planet forming disks 

The high metallicity of giant planets in our own solar system as well as planet for- 
mation models suggest that the atmosphere can be significantly polluted by evaporating 
planetesimals during the late stages of planet formation. The final composition of a plan- 
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etary envelope then depends not only on the radius at which the atmosphere is accreted, 
but also on the relative importance of gas and planetesimal accretion during atmosphere 
build-up. We parameterize these contributions to the relative abundance of an element X as 



"^X/H /x,solid Mjolid 



"X — stellar ~ ~r, m ^ \^ ~ /x,solidJ, 
'"'X/H ./s/g -'"gas 

SO that ax is the deviation of X=C or O from stellar. Here, m^J^ and rrf^^^^ arc the mass 
ratios of the element X with respect to H in the planetary envelope and star, respectively, 
/x,soiid is the fraction of X bound up in solids, /s/g is the mass ratio of solids versus gas in the 
disk, Mgoiid is the atmospheric mass accreted from solids, and Mgas is the atmospheric mass 
accreted from gas. In the presented model /g/g is set to 0.01, the observed grain to gas ratio 
in the ISM. The equation assumes that the amount of H bound up in solids is negligible and 
that the solids and gas are accreted from the same disk region, i.e. it ignores both planet 
and planetesimal movement. 

The effect of solids evaporating in the atmosphere, whether from accreting icy bodies 
or from core dredging, on the C/0 ratio is illustrated in Fig. 3 for an atmosphere accreted 
between the CO2 and CO snow lines (/c,soiid = 0.35 and /o,soiid = 0.67 from Table 1). Super- 
stellar C/0 ratios are achieved when Mg/Mg < /g/g, independent of the absolute value of /g/g. 
In Fig. 3, super-stellar C/0 ratios are produced when less than 1% of the total atmosphere 
mass comes from evaporation of solid material since /g/g = 0.01. Our parameterization also 
results in a prediction of the C/H ratio as a function of gas/solid atmosphere accretion. The 
C/H ratio is sub-stellar when the atmosphere is dominated by gas accretion (>99%) since 
the amount of H bound up in grains and grain mantles is neghgible compared to the gas H2 
content, while large amounts of C and O are in the grain mantles at these disk radii. 



3. Model extensions: disk and planet growth dynamics 

So far we have assumed a static disk, in which the gas and solid compositions are set 
purely by the evaporation lines of different grain mantle molecules, and where the entire 
planetary envelope forms at the same location. These assumptions ignore several dynamic 
processes active during planet formation. First, growing planets likely move within their 
natal disks during the formation process. While run-away gas accretion is simulated to be 
fast enough to effectively take place at a single radius, planetesimals may accrete over longer 
times. The solids that build up the atmosphere through either planetesimal accretion or core 
dredging may therefore originate at different locations compared to where the gas envelope 
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is accreted. 



Second, small solids drift inwards while planets are forming. After crossing a snowline, 
planetesimals evaporate over a finite length scale, which depends on a combination of the 
surface evaporation rate and the drift rate, enhancing th e gas phase abundance in the evap- 
orating species close to the snowline (IBrauer et al.l 120081 ). The C/H ratio incorporated into 
the atmosphere of a planet forming just interior of the snowline may therefore be higher 
compared to in a stationary disk, even super-stellar. Figure 3 includes predictions for the 
C/H and C/0 ratios for a planet that is forming in a region where the CO/H2 ratio is 
enhanced by a factor of two due to this kind of evaporation front. 

Finally the disk temperature profile is not static, but rather deperids on both the evolving 
lumin osity of the accreting star and the disk opacity (ILecar et al.ll2006l : iMakalkin fc Dorofeeva 
2OO9I ) . Therefore, even if a planet remains stationary it may accrete its envelope from gas 
and solids with evolving compositions. As discussed in §2 this should be a minor effect 
during gas accretion, but it may affect the planetesimal accretion pr ocess, which occurs over 



longer timescales. Fina lly, ice and gas also evolves chemically (e.g. iDodson-Robinson et al. 



2003 : IVisser et al.ll201ll ) during the disk life time. This is likely to only have a minor effect 
on the composition of forming planetary envelopes, however, since the abundances of the 
most important molecules H2O and CO are predicted to stay the same on 10® year timescales 
jAikawa et al.lllQQol : IVisser et allboill ). 



Ideally all these dynamic effects should be included when simulating the build-up of a 
planetary envelope quantitatively. Without these processes included in the model, we can, 
however, still predict some atmospheric composition trends based on how and where a planet 
envelope forms with respect to different snowlines, because of the short timescale of runaway 
gas accretion compared to most other disk processes, including migration, temperature vari- 
ations and chemical evolution. 



4. Discussion 

The model prescription above shows that protoplanetary disk models generally contain zones 
in the planet forming region where the gas is enhanced in carbon with respect to oxygen. 
The resulting planetary envelope composition will depend on this gas C/0 ratio as well as 
the amount of icy solids incorporated into the atmosphere through planetesimal accretion 
or core dredging, and the disk location where these solids form. As observations and our 
understanding of the relationship between atmosphere spectra and the overall atmosphere 
composition improve, the combination of C/0 and C/H may therefore provide an important 



- 7- 



probe of a planet's formation history. 

In summary, stellar C/0 is expected for planets forming through gravitational insta- 
bilities where all material is mixed (reduced C/0 and enhanced C/H are possible from late 
planetesimal accretion), and for planets forming interior to both the water snow-line and the 
carbon-grain evaporation line, where all volatiles are in the gas phase. Sub-stellar or stellar 
C/0 and super-stellar C/H indicates that a large amount of icy solids 'polluted' the atmo- 
sphere following gas accretion. Superstellar C/0 and C/H may be due either to gas accretion 
close to the CO or CO2 snowlines, or indicate that the atmosphere was substantially pol- 
luted by carbon-grains that have been postulated to be abundant between the carbon-grain 
sublimation line and the H2O snow-line. Finally, superstellar C/0 and substellar C/H are 
only consistent with atmosphere formation from mainly gas accretion outside of the water 
snowline. We note that this last case constrains the location of atmospheric accretion, while 
for super-stellar C/H, we constrain the formation location of the dominant planetesimal pol- 
lutants with respect to the locations of different snowlines in a disk. As both ice abundances 
and exo-planetary abundances of N-bearing species become better determined, C/N ratios 
will provide additional constraints on the planetary envelope formation history due to the 
complementary evaporation histories of C and N-bearing molecules. 

Testing these scenarios relies on well- determined atmospheric compositions, which are 
difficult to achieve even in our own solar system; the measured C/0 value for Jupiter is 
high, but this is believed to reflect a layere d atmosphere rath er than a constraint on the 



bulk composition of the planet's envelope (lAtreya et al.ll2002l ). since in a cool gas giant 



such as Jupiter, much of the 0-rich v olatile material originally present in the envelope may 



have condensed (IVisscher et al.ll2010l ). The super-stellar metallicities of Jupiter and Saturn 



are better established, an d indicate that larg e amounts of icy solids have been incorporated 



into their envelopes (e.g. iMousis et al.ll2009l ). An overall high C/0 ratio is therefore not 



expected, since accretion of solids outside of the H2O snowline decreases the C/0 ratio 
(Fig. 3). Current observations cannot constrain the bulk C/0 ratios of Kuiper belt objects 
(KBOs). However, our model predicts that KBOs formed beyond ~40 AU will have solar 
C/0 ratios, while those formed interior to ~40 AU and transported outward will have sub- 
solar C/0 ratios. We note that moons may not have formed at the same temperature as 
their host planets, complicating their use as chemical tracers. 

Because of the poor mixing of volatiles in gas giants in our own Solar System, th e 



compositions of hot Jupiters, where circulation is more efficient (IMadhusudhan et al.ll201ll ). 
may be comparatively easier to constrain, despite the technical difficulties involved. Of the 
handful of spectroscopically observed hot Jupiters, WASP 12b and possibly HD 189733B 
are candidates to have atmospheres originating from mainly gas accretion beyond both the 
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H2O and CO2 snowlines (iMadhusudhan et al.ll201ll : ISwain et al.ll2009l . I2OIOI ). pending con- 
firmation of the reported C/0 ratios and C/H ratios and tlieoretical confirmation of tfiat 
tfie observed compositions are representative of tfie bulk envelope. In both cases the planets 
must have migrated from their gas accretion zones beyond the snow-line into their current 
sub-AU orbits. 

In contras t, the directly ima ged young planets observed at 24-68 AU around the A 
star HR 8799 (IMarois et al.l 120081 ) may have formed close to their present location. Some 
of the planets are currently located in the exact region where the gas phase C/0 ratio is 
expected to the largest around an A star (Fig. 2). A discovered high C/0 ratio would 
be evidence for that the planetary envelopes formed through mainly gas accretion close to 
their current location. In contrast, stellar C/0 ratios in the planetary envelopes could be the 
product of either gravitational collapse or gas and planetesimal accretion further from the star 
followed by inward migration. Sub-stellar C/0 in the envelope would instead be evidence for 
either gravitational collapse or core accretion, followed by significant planetesimal enrichment 
interior to the CO snowline (Fig. 3). Spectra of HR 87 99b, situated at 68 AU, show a lack 
of CH4 absorption, possibly indicating a low C/0 ratio ( iBowler et al.ll2010[ ). Information on 
CO and H2O are however needed to co nstrain the atmosp heric composition further. With 
effective temperatures of ~800-1100 K (IMarois et al.ll2008l ). the young planets orbiting HR 
8799 are hotter than Jupiter but cooler than hot Jupiters. Whether spectra for these objects 
constrain their bulk envelope compositions is uncertain and merits further investigation. 

In general, accretion of oxygen depleted gas is a simple way of forming planet atmo- 
spheres with high C/0 ratios that directly arises from the current planet formation paradigm, 
consisting of core accretion in a disk with a thermal gradient. To predict atmospheric spectra 
based on the atmosphere formation mechanism requires, however, both detailed simulations 
that take into account gas and planetesimal accretion as the planet core moves through the 
disk, and further exploration of the circumstances under which bulk envelope compositions 
can be constrained from observable planetary atmospheres. 
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Fig. 1. — The C/0 ratio in the gas and in grains, assuming the temperature structure of a 
'typical' protoplanetary disk around a solar-type star (Tq is 200 K, and q — 0.62). The H2O, 
CO2 and CO snow-hnes are marked for reference. 
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Fig. 2. — The predicted gas phase C/0 ratio as a function of radius for five representative 
disks, ordered by spectral type, compared with the 'typical' disk model in Fig. 1. The 
derived temperature profile parameters, Tq and q, are listed. The C/0 ratios are calculated 
assuming that the stellar C/0 ratio is solar, i.e. 0.54, and a static disk. 
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Fig. 3. — The predicted atmospheric deviations of C/H (ac, red) and C/0 (ac/ao, black) 
from the stellar values, as a function of the fraction of the planet atmosphere mass that comes 
from planetesimal accretion (as opposed to gas accretion) for a planet forming between the 
CO2 and CO snowlines. The dashed lines assumed CO/H2 gas enrichment by a factor of 
2 because of an evaporation front close to the CO snowline, while the full lines assume the 
equilibrium CO/H2 gas ratio in a stationary disk. The stationary disk abundance ratios 
coincide with stellar abundances at 0.01 because this is the assumed grain/gas mass ratio 
fs/g in the disk. 
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Table 1: Evaporation temperatures and abundances of O and C in different forms witli 



rogen. Adopted 


model values in 


parentheses 




Species 


evap 


no 


nc 




(K) 


(IO-^xuh) 


(IO-^xuh) 


CO 


18-22 (20) 


0.9-2'^ (1.5) 


0.9-2^ (1.5) 


CO2 


42-52 (47) 


0.6^ 


0.3^ 


H2O 


120-150 (135) 


0.9^ 




Carbon grains 


>150 (500) 




0.6-1.2= (0.6) 


Silicate 


-1500 (1500) 


1.4= 





^The range of temperatures for ices corresponds to gas densities 10* -10^^ cm 

disk midplanes. 

^From ice and gas observations toward the CBRR 2422. 8-3423 disk 



suitable for 



^The range corresponds to estimates of organic content (jPraind 120031). The lower value is 



PontoppidanI l2006l ) . 



adopted to obtain a solar C/0 ratio. Silicate abundance is 1.2 from IWhittetl (I2OIOI ) and 1.4 
takes into account the additional refractory component. 



